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NONLINEAR STOCHASTIC PARTIAL DIFFERENTIAL 
EQUATIONS WITH SINGULAR DIFFUSIVITY AND GRADIENT 

STRATONOVIGH NOISE 

lOANA CIOTIR AND JONAS M. TOLLE 


Abstract. We study existence and uniqueness of a variational solution in 
terms of stochastic variational inequalities (SVI) to stochastic nonlinear dif¬ 
fusion equations with a highly singular diffusivity term and multiplicative 
Stratonovich gradient-type noise. We derive a commutator relation for the 
unbounded noise coefficients in terms of a geometric Killing vector condition. 
The drift term is given by the total variation flow, respectively, by a singular 
p-Laplace-type operator. We impose nonlinear zero Neumann boundary con¬ 
ditions and precisely investigate their connection with the coefficient fields of 
the noise. This solves an open problem posed in [Barbu, Brzezniak, Hausen- 
bias, Tubaro; Stoch. Proc. AppL, 123 (2013)] and [Barbu, Rockner; J. Eur. 
Math. Soc., 17 (2015)]. 


1. Introduction 


We consider existence and uniqueness of solutions to the following (multi-valued) 
nonlinear Stratonovich stochastic diffusion equation in 


( 1 . 1 ) 




N 

dXt e div [sgn (VAj)] dt + o dPl, 

i=l 

Xo = X, 



in (0,T) X O, 
inO, 

on (0,r) X do, 


where O is an open, bounded domain in d > 2, with (sufficiently) smooth 
boundary such that O or dO is convex. Here, for A > 1, 6^ : 0 — >■ K'^, 1 < * < A 
are “coefficient fields” and /3 = (/3^,... ,P^) denotes an A-dimensional Brownian 
motion on a filtered (normal) probability space (h, A, , p) . The initial 

datum is chosen as a: G L‘^{0), or, more generally, as a; G L^(H, Jq, P; L^((!l)). 
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Here, v denotes the outer unit normal on dO. The multi-valued graph ^ i-)- sgn(^) 
from into 2®“^ is defined by 


sgn(0 := 


i^r 

{CeK^IICI <1} 


, if e ^ 0, 
, ife = o, 


for all ^ Because of the multi-valued diffusivity term, the equation becomes 

formally a stochastic evolution inclusion, as have been studied e.g. in [M1I32II38] . We 
denote by I'l the Euclidean norm of R'^, and by (•, •) the Euclidean scalar product 
of R'^. 

Set 


( 1 . 2 ) 



: O 


and denote by b* its transpose. We have that equation (HID is formally equivalent 
to the Ito stochastic partial differential equation, 

(1.3) 

dXt £ div [sgn (VXj)] dt+^ div[b*bVXt] dt + (bVXj, djdt), in (0, T) x O, 

< Xq = X, in O, 

-^ = 0. on(0,T)xdO. 

on 

A similar equation was studied in [7] for the case of a dissipative drift, using the 
method of Brezis-Ekeland’s variational principl43- On the other hand, equations 
with singular drift of the same form have been studied in [S UIOll^ for additive and 
multiplicative bounded noise, respectively. See |37j for a multiplicative Stratonovich 
stochastic equation with a similar drift term. Those results do not apply to our 
case since the noise coefficient 


(1.4) u !->■ (bVu, •) 

is not bounded on the state space {O). In [3D], existence and uniqueness as well 
as regularity have been investigated for the stochastic mean curvature flow with 
unbounded noise. The methods used are related to ours, however, the structure of 
the equation prevents a direct application to our situation. 

Additionally, in our main Theorem 14.11 we will derive existence and uniqueness 
results also for the singular p-Laplace equations with p £ (1, 2), 


N 


dXt = div [\XXt\P-^VXt] dt + VVt) o dpi, 

in (0,r) X O, 

Xq = X, 

in O, 

dXt 

on (0,T) X do 

dn =0- 


^Equation II1.3II with Dirichlet boundary conditions (instead of Neumann boundary conditions) 
is also being investigated in m- However, the preprint of m became publicly available after our 
revised work was submitted for publication. We point out that the method used in m is different 
from ours. 
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Due to the lack of strong coercivity of the drift operator, we shall employ so- 
called stochastic variational inequalities (SVI), with the aim to construct solutions 
to (HHl) in a weak variational sense. Even for bounded noise, singular equations of 
the above type are generally not known to satisfy an ltd integral equation — not 
even in the (analytically) weak sense. Compare with [8l|9l|30l[33] for related works 
employing SVI-frameworks. Using a rough path approach, equations with similar 
noise were studied in |171I29] . A similar equation with linear drift is investigated 
in m- We would like to point out, that the solutions of the work at hand are strong 
solutions in the probabilistic sense, meaning, in particular, that the solutions are 
functions of the given Brownian motion. 

The natural energy space for the (Neumann) total variation flow, the p-Laplace, 
respectively, would be BV{0), the space of bounded variation functions, respec¬ 
tively, the Sobolev space W^’^{0). However, on the level of approximations, we 
shall work on the smaller space H^{0). One reason is, that we are using viscosity 
approximations, namely, we are adding a regularization term eA, and taking e 0. 
In particular, this allows us to consider the gradient-type SPDE for the borderline 
case of a monotone drift operator (p = 1) which cannot be treated within the scope 
of reflexive Gelfand triples {p > 1), as e.g. has been done in [7l[T0] for Dirichlet 
boundary conditions. 

Another property, necessary for our arguments, is the mutual commutation be¬ 
havior of the diffusion coefficients, as well as the question of commutation with the 
Neumann Laplace — in order to obtain these, we introduce a condition from dif¬ 
ferential geometry, similar to the notion of Killing vector fields, see AssumDtion l2.ll 
and Appendix lAl below. In this context, we prove that, under our assumptions, the 
first-order partial differential operator dm), which corresponds to an infinitesimal 
vector held action, preserves Neumann boundary conditions, see Lemma l2. 81 below. 

According to [46], the interest in studying this type of equation comes from its 
use for simulations in the tomographic reconstruction problem, which has several 
applications, for instance in medical imaging and general image processing. 

More precisely, the binary tomography methods are proposed in [31] as a simpler 
inverse problem of reconstruction. Being still an ill-posed problem, it needs to be 
regularized, and this may be done for instance with the total variation (T.V.) 
regularization. In order to numerically solve the problem, a fast and efficient 
T.V./L^ minimization algorithm based on the “Alternate Direction of Minimization 
Method” (A.D.M.M.) has been proposed in [T1[S2]. Finally, a singular stochastic 
diffusion equation with gradient dependent noise is used to rehne the solution ob¬ 
tained by the A.D.M.M. algorithm, see also the related Example 12.31 below. The 
time dependent (deterministic) T.V. image restoration problem has been studied 
e.g. in [15]. We refer to [34] and the references therein for a stationary stochastic 
approach. 

Therefore, the present work gives rigorous theory to support the use of this kind 
of equation for numerical results such that those in [46] . However, the authors 
of [46] are posing the problem for an Ito-equation instead of a Stratonovich one, 
see also m- 

Another possible interest of studying stochastic differential equations perturbed 
by this type of noise comes from the applications in modes of turbulence (see m)- 

Discussion of an approach via transformation. Following the classical works 
[26l|48] , we can also think of an alternative access to our equation, which, however. 
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must fail even on a heuristic level. Here, we shall briefly discuss this approach and 
point out the difficulties. 

Let y € LF'{0) and consider the following deterministic PDE 


dYt{0 G div[sgn(Vrt(e))] dt, Yo{0 = te{0,T],^eO, 

where we impose Neumann boundary conditions. For initial datum y G 
a unique weak solution in the Gelfand triple C C H~^ was constructed 
in [32l Theorem 2.6]. For initial conditions y G L'^{0), see [21131123]. 

Let b be as in (11.21) . and assume merely that bi € C^(0; for 1 < i < TV. For 
t G 10,T], ^ G O, uj G fi, define 

(1.6) := yt(^ + MO = y(0- 

A similar transformation approach can be found in m for linear equations and 
in [5^ for the case of conservation laws. See also [T3] for other nonlinear SPDFs 
treated by this transformation. 

Assume for a while, that we have a pathwise Ito formula available (that is, for 
uj G ft, fixed), ignoring the lack of regularity of {x,t) Yt{^ + b*(^)a;) =: F{x,t) 
for a moment: 


F{Pt,t) = F{0,0) + 


N 

E 


d,,FiPs,s) odPl 


f 


dtF{/3s:S)ds, 


see [nilll|27l|28]. By the chain rule, we would obtain that for dt-a.e. t G [0,r], 
possibly outside an exceptional subset of O, 


N d 

Xt&y+ div[sgn(VXs)]fis +VV / b^d^^Xs o d/dl, 

Jo i=i Jo 

which is a pathwise representation of equation (HI). The Stratonovich correction 
term is formally given by 


d N 


^[vm-),/3]* = ^EE / VXs))(is = i f div[b*bVA:s 

1 __1*^0 'J 0 


ds, 


k—1i—1 


where t i-A [•, •]* denotes the quadratic covariation process, compare with [271128] . 

Fven if one finds a way to deal with the measurability issues, the direct ap¬ 
plication of this approach must fail due to the lack of regularity, since, according 
to [T8l[T9l[25] , good Holder estimates for the solution (and for the gradient of the 
solution) to the parabolic p-Laplace equation usually hold only if p > thus 
sorting out the total variation flow. 


Organization of the paper. After a brief part on notational conventions of this 
work, we shall give our assumptions and discuss the resulting properties of the 
noise coefficient operators m in Section [2]— in particular, we establish the com¬ 
mutation relations which we shall need subsequently. Our notion of SVI-solutions 
(to equations with gradient-type multiplicative Stratonovich noise) is provided in 
Section 131 In Section 0] we shall first derive a useful a priori estimate in 
and after that go through several approximation steps necessary for proving the 
existence of a solution. The uniqueness of SVI solutions is proved in Subsection 
14.21 For the reader’s convenience, we shall provide some results on Killing vector 
fields in the appendix. 
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Notation. We shall recall a few standard definitions and fix notation which will 
be used later. 

We set V := H^{0) = the standard first order square integrable 

Sobolev space and H := the Hilbert space of (classes of) square integrable 

functions with respect to the Lebesgue measure. We also consider the second order 
square integrable Sobolev space H‘^{0) = W‘^'‘^{0). We shall write 
and so on, if the context is clear. We denote the inner product in H by (•, ■)h- V* 
denotes the topological dual of V with dualization denoted by (•, •). Let 
be the usual first order p-integrable Sobolev space. For u G L^{0) we define the 
total variation semi-norm by 

IIuIItv := sup a u div 77 

and let BV be the space of functions of bounded variation, that is, 

BV[0) := {u G L\0) I ||u||tv < 00 }. 

For a proper, convex, lower semi-continuous (l.s.c.) function ^ : H ^ [0,-l-oo], 
we denote the subdifferential by 9$. The graph of 9$ consists precisely of the 
pairs of elements {x,y) £ C H x H that satisfy {y,z — x)h < for 

all z £ H. In this context, we may also write y £ 9$(re), where we identify the 
subdifferential as a multi-valued map 9$ . H 2^. 

We say that a function X £ L^([0,T] x is {J’tj-progressively measurable 

if is H([0,t]) ( 8 ) J^t-measurable for all t £ [0,T]. The domain of (unbounded) 

linear operators A is denoted by dom(H), and by the same notation, we denote the 
effective domain of convex functionals or multi-valued graphs. By C, we denote a 
positive constant that may change its value from line to line. 

2. Hypotheses and commutation relation 

Suppose that O C is a sufficiently smooth, open, bounded domain. Denote 
the surface element on dO by Denote by u the outer unit normal on dO. 

Below, we collect our assumptions on the “diffusion matrix” b and prove some 
essential properties of the associated partial differential operators. Briefly summa¬ 
rized, we are assuming conditions to ensure that 

• the first-order partial differential operators associated to the rows of b are 
well-defined unbounded skew-symmetric linear operators on L^{0), see 
Assumption 12.11 (i) and Lemma 12.61 below; 

• the groups of diffeomorphisms generated by the rows of b mutually com¬ 
mute, see Assumption 12.11 (ii) and Lemma 12.71 below; 

• the partial differential operators associated to the rows of b leave the 
domain of the Neumann Laplace invariant and commute with its resolvent, 
see Assumption [2T] (i), (hi), (iv) and Lemmas 12.81 and 12.101 below. 

We note that the commutation assumptions are typical for gradient-type noise, 
even for linear stochastic equations, see and dU Section 6.5]. 

Assumption 2.1. Suppose that the diffusion coefficients bi £ (^^(O;R'^), 1 < j < 
N, satisfy the following conditions: 

(i) (bi, v) = 0 on do for all 1 < i < N. 

(ii) Either N = 1, or bfdklrl = b^dkb{ on O for all 1 < k,j Si d, 1 < i, I < N, 
i ffl. 


y£C^{o-X)i hlUo 


(C>;R<i) ^ 1 
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(iii) div6i = 0 and {Abi,bi) — 0 on O for all 1 < i < N (where the Laplace 
operator acts componentwise). 

(iv) {(ybi,v),bi) + {{\Ibi,hi),v) = 0 on dO for all 1 < i < N (where the 
gradient acts componentwise). 

By Lemma IA.2I in the appendix, sufficiently smooth vector fields bi that satisfy 
(iii) and (iv) above, are precisely the so-called Killing vector fields, see (lA.ll) in the 
appendix for the definition. 

Remark 2.2. Condition (i) and (iii) in the above Assumption sort out any examples 
with stochastic perturbation for the case d = 1. Indeed, let O = / be a bounded 
interval, so that clearly ly = ±1 at the endpoints of I. One the one hand, condition 
(i) implies that 6 = 0 on dl. On the other hand, div b — b' = 0 implies that b must 
be constant on I. Hence 6 = 0. 

Altogether, condition (i) ensures that the noise coefficients respect Neumann 
boundary conditiont0, see Lemma [2.81 below. 

Example 2.3. Let fV = 1, d = 2. Let O = {C e | |CI < R}, R > 0- Let 
b{f) := ($ 2 ,—Cl)- Then (11.11) becomes 

t dXt € div [sgn (VAt)] dt + (^ 2^1 At - £,id 2 Xt) o dft, in (0, T) x O, 

I Aq = X, in O, 

1 dJC 

= 0 , on (0,r) X do. 

^ ah' 

Example 2.4. Let A = 1, d = 3. Let O = {C G |Cl < R}, R > 0- Let 

KO ■= (6 -6,Cl -C 3 ,C 2 - Cl) and denote 1 := (1,1,1) (clearly, 6(C) = C x 1 ). 

Then 6 is a Killing vector field and dnj becomes 

f dXt G div [sgn (VAt)] df + (C X VAt, 1) o dft, in (0, T) x O, 

I Aq = X, in O, 

1 

^=0, oii(0,T)xdO. 

^ ou 

One can replace 1 by any constant vector Co G \ {0} and get that 6(C) := C x Co 
still satisfies Assumption 12.II 

Remark 2.5. Note that: 

(i) The above vector fields C (C 2 ,—Ci) and C C x Co resp. are the 

infinitesimal generators of the rotation groups SO{2) and SO{3) resp., see 
e.g. [35]. They generate groups of rotations around the origin, leaving 
balls centered at the origin invariant, which explains why the respective 
domains are chosen as above (Co spans the axis of rotation). 

(ii) Let d = N. The example of constant vector fields bj = Sij are precisely 
the infinitesimal generators of groups of translations (violating Assumption 
12.11 (i) on balls). The d-torus leaves the translation groups invariant, 
and is still a bounded, convex domain, leading either to periodic boundary 
conditions or to a setting for compact manifolds without boundary. 


^If b=l on 0 = 1 = {0,2w), then ^ i—cos^ has Neumann boundary conditions, however, 
b- (cos 5)' does not. 
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Recall that the domain dom(—A) of the Neumann Laplace in the weak sense is 
given by all elements u S H^{0) such that Au G Lp'{0) and such that 

f vAud^ = — f (Vu,Vit)d^ 'ivGH^{0). 

Jo Jo 

For u G dom(—A), the normal derivative ^ belongs to (being the dual 

of the space of traces H^^‘^{dO)) and is zero, see e.g. p. 250] for details. As we 
assume smooth boundary, the normal derivative is given by ^ = (Vu, v) S'‘^“^-a.e., 
whenever u G C'^{0). Hence, 

(2.1) C:={uG C^(0) I (Vu, i/) = 0 S^'-^-a.e.} 

is a core for the Neumann Laplace, that is, dense in dom(—A) w.r.t. to the graph 
norm 

\MLmi-A) ■■= [ (lAup + luHde 

Jo 

On the domain H^{0), we define the linear operators Bi, 1 < i < N as 

B, : H\0) L^{0) 


B,{um ■■= 

= diY[hiOum, VueiLi(O), 

where bi satisfies Assumption 12.11 

Lemma 2.6. Assume Assumvtion \2. 1\ Let us collect the following properties: 

(i) The space H^{0) is the domain of skew-adjointness of Bi, 1 <i < N, that 


IS 


B,u = -B*u, l<i<N, uGH^{0), 

where B* denotes the adjoint operator in L^lfD). 

(ii) For all u G H^(0), it holds that 


( 2 . 2 ) 


L 


uBiU df = Q, 1 < i < N. 


Proof. Let u G H^{0) and fix 1 < z < iV. 

(i) : By the Gauss-Green theorem, for v G H^{0), taking Assumption 12.11 fil 

into account, 

{B*U,v)l:‘{0) = {u,B,v)l2(^0) 

= / u{bi,Vv)d^ 

Jo 

= — diY{biu)v d^ - 

Jo 

= - {bi,Vu)vd^ 

Jo 

= -{B,U,v)L-i( 0 )- 

The density of H^{0) C Lf{0) yields (i). 

(ii) : This follows directly from (i). 


uv{bi, u) dS 


d-l 


IdO 


□ 
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For 1 < * < fixed, let L'^{0), t G M., denote the Co-group of 

linear operators associated to Bi, such that, in particular, 

= BiU, ueH\0). 

t=o 

Lemma 2.7. The groups , 1 < i < N, t G R mutually commute, whenever 
Assumvtion \2. 1\ holds. 



Proof. For TV = 1, there is nothing to prove. Let bi, bi, i ^ I be as above. Let 
/ G C^{0). Define the commutator [bi, bi]f := BiBif — BiBif. By Leibniz’s rule, 

[h,bi]f= Y. b':d,ibid,f)-bfdkmf) 

l<k,j<d 

= Y b'^bjdkd.f + b'ydkbid,f -bfbldkd.f -bfdkbld,/ 

l<kj<d 

= Y {b'^dkbi-bfdkbDd.f 

l<k,j<d 

= 0 , 

where we have used Assumption 12 .1 1 fii'l in the last step. 

Now, for any 1 < z < N, denote hy Zl : O ^ O, t G [0, oo) the flow of 
diffeomorphisms on O corresponding to the vector field action of bi, that is, 

j^Zl = b,{Zl), t>0, Zl=^GO. 

Then, since we have proved above that \bi,bi] = 0 for any i ^ I, and from the 
property of the vector fields to be divergence free, we get that the Z"^, 1 <i < N 
mutually commute (in the sense of composition of maps), see [331 Ch. 1.2, Exercise 
3]. However, it holds that (e‘^‘u)(^) = u{Zl{£^)) and one easily deduces that the 
groups of operators commute on L^. □ 


Lemma 2.8. Suppose that Assumption {ZJi holds. Let u G C^{0) be a scalar 
function with (Vu, p) — 0 on dO. Then it holds that (V((6i, Vm)), z/) = 0 on dO 
for every \ < i < N. 


Proof. Fix 1 < z < and set b := bi. Let zy G C^{0) be a test-function. We claim 
that for any u G C^(0) with (Vu,v) = 0 on dO, it holds that 

(2.3) [ ry(V((6, Vzz)), p) dS^^-^ =0 Vzy e C^{0). 

JdO 

In order to prove (12.3p . we first apply Gauss’s divergence theorem to the vector 
field F := r]y{{b, Vzz)), and get that 


{r]V{{b, Vu)),v) dS' 


d-l 


div Fd^. 


JdO Jo 

However, div F — (Vzy, V((6, Vzt))) -I- riA{{b, Vzz)). Let us begin with investigating 
the second term. By the Killing assumption, we have the commutation on suffi¬ 
ciently smooth functions (cf. Theorem IA.3l in the appendix), thus, riA{{b,\/u)) = 
r]{b,S/Au) = {rib,\7Au). Integrating by parts, we get that 


f {rib,\7Au) d^ = — f div{r]b)Aud^ + f 'q{b,v)AudS'^ 

Jo Jo JdO 
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The latter term is zero by {b, v) — 0. Also, since both 77 and b are smooth up to 
the boundary, div(? 76 ) G H^(0). So we can use the Neumann boundary condition 
for u to get that 

— f div{r]b)Aud^ = ( (V div( 77 &), Vm) 

Jo Jo 

Clearly, as div 6 = 0 , we get that V div( 77 &) = V({&, V??)). Hence 

[ div FdC= [ [(Vr7,V((&,VM))) + (V({&,V77)),VM)]de, 

Jo Jo 

differentiating out this term yields 

[ div Fd^= [ [{{Db)-Vu,Vf]) + {{D^u)-b,Vr]) + {{Db)-Vr],\/u) + {{D^r])-b,S/u)]d^, 

Jo Jo 

where, Db denotes the Jacobian of b and denotes the Hessian of a scalar function, 
denotes matrix multiplication. 

However, Db is skew-symmetric with respect to the Euclidean scalar product 
due to the Killing assumption, see (lA.ll) in the appendix. Hence 

{{Db) ■ Vm, Vt?) = -{{Db) ■ V??, Vm) 

and the above term becomes 

[ div Fd^= [ [{{D^u)-b,\/r]) + {{D^'n)-b,Vu)]dC 

Jo Jo 

With Einstein’s summation convention, interchanging the order of differentiation, 
f [didjuDdiT] + didjrjDdiu] d^ = f [djdiuDdirj + didjr]Ddiu] d^. 

Jo Jo 

Integrating by parts in the first term yields 

f djdiuDdiTjd^ = — f diudj{DdiT]) d^ + f diuDdi’qvj dS'^~^. 

Jo Jo JdO 

Now, in the boundary integral term, we can separate the sums over j and i resp. 
and get that this term becomes zero by (5, v) = 0. Furthermore, 

—diudj{D diTj) = —diudjWdiT] — diubJdjdir] = —diuVdjdirj, 
as we have that div 6 = 0. Finally, the remaining terms cancel, and we get that 


0= [ div FdC= [ 77(V((6,Vu)),p)dS''^-i \fr]GC^{0). 

Jo JdO 


□ 


We shall need the following commutation result. Denote the resolvent of the 
Neumann Laplace by Js '■= (Id —<5A)~^, J > 0. 

Theorem 2.9 (Shigekawa). Fix 1 < i < N . Suppose that there exists a linear 
subspace D C dom(—A) such that the following conditions hold: 

(i) A{V) C dom(Hi), 

(ii) Bi{V) C dom(-A), 

(hi) D is a core (see for the terminology) for (—A, dom(—A)), 

(iv) dom(—A) C dom(i?i) and dom(—A) C dom(H*), 

(v) For any u GD, it holds that 

BiAu = ABiU. 
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Then for all 5 > 0, and every u G doia(Bi), it holds that 


Bi — J§BiU. 

Proof. See [JU Theorem 3.1 and Proposition 3.2]. □ 

Lemma 2.10. Assumvtion \2. 1\ implies that all of the conditions of Theorem \2.t]\ 
are satisfied for all 1 < i < N. 

Proof. Assume the conditions of Assumption 12.11 Fix 1 < i < A^. First note that 
Lemma [221 implies that dom(i?i) = dom(i3*) = H^{0). Let V := C°°{0) DC, 
where C := {u G C'^{0) \ (Vu^iy) = 0 S^^^^-a.e.}. Obviously, 21 is a core for 
dom(—A). Hence (i), (iii), (iv) are clearly satisfied, (ii) follows from Lemma [2.81 
The commutation on smooth functions (v) follows from Theorem lA.31 since bi is a 
Killing field by Assumption 12.11 and Lemma IA.2I □ 

Let us also define Bf : P[^{0) —> (22^(0))* by 

B^u:= -B*BiU, uGH\ 0), l<i<N. 

In the sense of Schwartz distributions, it holds that 

N 

'^B'^u = div[h*hVu]. 

Set S := H^{0). We thus have a Gelfand triple 

ScHcS*. 


3. Stochastic variational inequalities (SVI) 

Let P = (/3^,..., P^) be a TV-dimensional Brownian motion on a filtered (normal) 
probability space (with expected value E[y] = f^YdP, Y G 

L^(fl)). We consider the following SPDE on H = L'^{0), where O C is a 

smooth, open, bounded domain such that O or dO is convex, 

(3.1) 

, N N 


dXt G div [dt (VW)] dt+-J2 BjXt dt + J2 B^Xt dPl , in (0, T) x O, 

i=l i=l 

in O, 

= 0, on{0,T)xdO, 


Xo = X, 

dXt 
dv 


here, ip := dip C x is the subdifferential of 93 := i| • |^ for p G [1,2), which 
is multi-valued for p = 1, i.e. d(f i-A |C|)(-) = sgn(-). More precisely, after fixing 
p G [1,2), let 

;/o‘^(^“(0) AuGH\ 

-1-00, if M G \ H^. 


I>(m) := 


$ is a proper convex functional on but might fail to be lower semi-continuous. 
Let us define 

$(u) := cll>(u) := inf < liminf $(u„) | Un ^ u G Lp {O) strongly i , 

L n—¥oo ' J 
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^{u) := 


the so-called lower semi-continuous envelope of cf. [5j Proposition 11.1.1]. The 
l.s.c. envelope is given by, for p £ (1, 2), 

.= I fo if« e W^:P(0) n L^(0), 

\+oo Hue L^{0)\W^'P{0), 

and for p = 1, 

^ ||w||ti/ if u g BViO) n L'^iO), 

^+oo iiueL^{0)\BV{0)., 

where we suppress the dependence on p in the notation. Obviously, $ is convex and 
it is easy to see that $ is lower semi-continuous on H. Moreover, $ is Gateaux- 
differentiable in u with derivative given by 

D^{u){v) = [ {r],\7v)di, 

Jo 

with r]{^) e '${Vu{^)) for a.e. ^ e O. In fact, $ coincides with the lower semi- 
continuous hull of $ on if, and we have for it G that 

{—div ?7 I T] e rj e di(Vit), d^-a.e.} C d^{u). 

However, the full characterization of 9$ (already in the space L^{0)) is involved. 
We shall omit its precise characterization and instead refer to [S]. 

Equation (I33D is then written in relaxed form as 

N N 


dXt e —d^(Xt) dt -h — y ] BfXt dt + y ] BiXt dPl, t e (0, T), 


(3.2) 


Xr, = X. 


Motivated by 012], let us define our notion of a solution to (13.21) . 

Definition 3.1. Let x e T > 0. An {Jil-progressively mea¬ 

surable process X e L^([0,T] x H;L^((!1)) is called an SVI-solution to (13.21) if 

(i) (Regularity) 

$(A) G L^i[0,T] X H). 

(ii) (Variational inequality) For every Z G L^([0,T] x Vt\H^{0)) such that 

there exist Zq G G G L‘^{[Q,T] x Vt-L\0)), {Ft}- 

progressively measurable, such that the following equality holds L^{0), 
that is. 


(3.3) 


(3.4) 


Zt — Zo 


/ G,ds + -V / Bfz,ds + y2 / B^Z.dpl, 
Jq ^ .-i .10 .i io 


P-a.s. for all t G [0,T], we have that the following variational inequality 
holds true 

ft 


^E\\Xt-Zt\\h^a)+^l ms)ds 




■E 


^{Zs) ds 


-E 


/ {Gs,Xs — Zs)l'^(o) ds, 
Jo 


(3.5) 
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for almost alH G [0,r]. 

Moreover, if X G L^(n; T]; L^(0))), we say that X is a (time-) continuous 

SVI solution to (lO) . 

Remark 3.2. Practically, the test-process Z needs to satisfy Z G i^([0, T] xfl; H‘^{0)), 
we shall provide in (14.51) below that a process Z of the form (13.41) in fact exists (see 
also (14.41) below). 


Inequality (13.5p is obtained by formally applying the Ito formula for the square 
of the H-noim to the process 


N N 

d{X -Z) = (-5$(X) -G)dt+^Y^ B({X - Z) dt+ '^B,{X - Z) dp\ 

^ i=i i=i 

taking expectation and using the subdifferential property. 


4. Existence and uniqueness 

4.1. Existence. 

Theorem 4.1. Let x G Xq, P; i?). Then there is a unique continuous SVI 

solution X G L'^(LI;C{[0,T]; H)) to (13.21) in the sense of Definition \3.1[ For two 
SVI solutions X, X with initial conditions x,y € L^(n, Xo, P; i?), resp., we have 

esssupE||Xt - YtWjj < E||a; - y\\jj. 
te[o,T] 

Proof. Recall the notation H = S = . We first assume an initial condition 

xq G L^(fI,Jo,P;S') and, in the last part of the proof, we shall generalize to xq G 

In order to prove the existence of the solution, we need to take a threefold 
approximation for equation dUD. Therefore, we consider the following regularized 
equation, 

(4.1) 

dXf = Js div 4-^ (v JaXf dt + eAXf’^’^dt 

^ ^ XI ^ / X l\ 

+x E(sf)2(xr^’')dt+E Bf dpi, 

) ^ i—l i—1 ^ ' 

Aq’ =xo, 

a 


where 4/^, A > 0, is the Yosida approximation of di (cf. [6l p. 37]), J^, 5 > 0, 
is the resolvent of the Neumann Laplacian L := —A, i.e., Js = (Id—(5A)~^ and 
Bf{-) = Bpjsi-)). 

By m Theorem 2.4], we have that the there exists a unique {Xt}-adapted so¬ 
lution with X G C'([0,T];iJ) D L^([0, T]; 5) P-a.s. such that (14.11) holds P-a.s. in 
L^([0, T]; S'*). We note that our hypotheses guarantee that the conditions needed 
for [3 Theorem 2.4] are satisfied. 

Step I (the estimate in (O)): 


in (0,T) X O, 

in O, 

on (0,T) X do, 
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Considering Jq,, ol > 0, the resolvent of the Neumann Laplace operator — A, we 
define the sequence of semi-inner products on H 

:= {{-A)aU,v)fj , u,v€H, 

where (—A)^ is the Yosida approximation of the operator —A, i.e., (—A)^ = 
i(Id—J q) = —AJa and the induced semi-norms 


R L := 


(-A)^u 


H 


u G H, 


where (—A)a denotes the operator square root. 

Since they are continuous on H and for all m G S' we have that 

ll«L—^ I|V^^IIl 2(0;R<*) asa^O. 

We shall apply the ltd formula [HI Theorem 4.2.5] to (14.1|) with the functional 
u i-A ||m||q, for A and 6 fixed, and we get that for all t G [0,T] and P-a.s. 


(4.2) 


Y, 


€,X,6 


ll^olla+2 [ ((-A)„YfY^,j,divvI/^(VJ 5 YjY^))^ds 

+2e [ s((-A)„X|’^-^AX|’^-^)s.ds 

+ 2 ^ / 


2=1 

N 


+ Y, / i-A)lBfXI’^’^ 


i=l ■ 


ds. 


H 


By well-known properties of the resolvent (as symmetry in commutation with 
the Yosida approximation) and keeping in mind that the the operators Bi commute 
with the resolvent of the Neumann Laplace by Theorem 12.91 and Lemma 12.101 we 
can easily see that, by setting, 

$ a ( w ) := [ ip^iVu)d^, uGS, 

JO 

and setting v = that 

((-A)ai;, Jidiv^''^ {\/Jsv))„ = --(p - JaV,d{^x o Js)v)h, 

"a 

cf. [IHl Proposition 11.7.8] for the chain rule for subdifferentials. By using the 
argument of [SS] Equation (3.7)jl (here, the convexity assumption on the boundary 
is needed, see also [32j Example 7.11], where the heat kernel estimates of 
are applied), we see that 

-{JaV - t,5($a o Js)v)h 
a 

<— i^\{JsJaV) — ^\{Jsv)) = — {^xiJaJsv) — ^\{Jsv)) < 0 . 
a a 


^For Dirichlet boundary conditions on piecewise convex domains, this result has been proved 
directly without the use of heat kernel estimates in in [9l Appendix]. 
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Note that, since v G H^{0), we have that 


s{{-A)aV,Av)s‘- < - ||(-A)a?;||^ 


To see this, just take into account that 

0 < - \7JaV,\/v - = {{-A)aV,AJaV)H -s{{-A)aV,Av)s>- 

Furthermore, by commutation (see Theorem 12.91 and Lemma r2.10L 

{{-AUv,{Bf)\)^ 

H 


f H 

= - ii-AUv, JsB:Bfv)^ = - {Bt{-A)^v,Btv)^ 

H 


(^Bf{-A)kv,Bf{-A)lv)^ = - ||(-A)|i?fn 


By going back and replacing in (14.21) we get that P ( 8 ) ds-a.s.. 


(4.3) 


a: 


€,X,S 


< \\xo\\l-2e [ ||(-A) 

Jo 


ds 


N 


ds 


/ (-A)2i?f(AjA^) 

i=i Jo 

+2^2 / {{-A)^XlAS,BfXlASdpi) 

.•_1 Jo 


H 


2=1 

N 


+ E/ (-A)I Bf 

.-1 Jo 


ds. 


H 


iaking the expectation 




(4.4) 


E 


va: 


e,A,(5 


L2(C9;If 


- 2eE 


AAg’ ’ llj:^ ds < E II Va;o||^2((p.]5 


Step II (5 -G OJ; 

We shall pass to the limit in (14.111 for d ^ 0 by using Theorem 2.2 from [7]. 
Note that Fatou’s lemma (after passing on to an a.e. convergent subsequence) and 
^xiJs ) < ‘bA(-) (which holds e.g. by [32l Example 7.11]) imply that $^ 0^5 —s- 
in Mosco sense as <5 0 (for the terminology, see [4]). Therefore, we have that 


Js div 4-^ {VJs (•)) + eA (•) A div 4-^ (V (•)) + eA (•), 

as (5 0 and so for the corresponding inverse subdifferential operator^ Also, it is 

clear that for v G S, 

X{Jsv) —)► Vu, strongly in L‘^{0;R‘^), 


as d —>■ 0, which is sufficient for the strong convergence of the Co-groups of linear 
operators associated to Bf to the Co-group associated to Bi, see e.g. [Ej. There¬ 
fore, we can apply [71 Theorems 2.2 and 2.3] (note that we do not need that the 


'^Note that, <I>" —>■ "I? in Mosco sense implies that (9<t> in G-sense and ((9<t>"') ^ —I (S'!?) ^ 

in G-sense, see [4l[6| . 
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semigroups converge in C^{[0,T]] H), as we do not assume any time dependence 
for our noise coefficients) and obtain that P-a.s. as <5 —)• 0, 

^e,A,5 ^ weakly in ([0, T]; S') and 

weakly*in L°° ([0,r];iJ). 

Combining with (I4.4L we get by weak lower semicontinuity of the norm that 


ess sup E 
ie[0,T] 


VX 


£,A 


L2(C);1 


■2eE 


AXf’^ll, ds < EllVxn 


S \\H — 


Il2(CI;1 


We have proved that there exists a strong solution (in the sense of 0 ) to 
dXt’^ = div 4-^ dt + eAXf ^ dt in (0, T) x O, 


(4.5) 


IN . N 

+ 2 E BfiXt'^)dt + E B,(X,^’^) dpl, 


•V'£, A 
^0 — ^ 0 ; 

“'Xo, 


i—1 


2=1 


dv 


in O, 

on (0,r) X do, 


for initial datum xq S S) which is of the particular form (13.41) as claimed in 

Remark 13.21 

Step III (X 0); 

By applying the ltd formula with it i-A i ||m||^ and the expectation to the differ¬ 


ence 


eAXf ^ - Gt) dt 


d(^Xt^-Zt) = (div(vl/^(vxf^)) 

N N 

+ 2 ^ dt + Y^B, {^XY - Zt) dPl 


2=1 


2=1 


for Z and G considered as in Definition 13.11 we see that X®A js ajgo a SVI solution 
to (14.ip . i.e. 


(4.6) 

Xl’^ - Zt 

+E [ [ (VXJA) d^ds 

H Jo Jo 

-l-eE 

r / (VXJ■^VXJ-^-VZ,)d^ds 

do JO 

— 2® Zo\\^ 

+ E f f ip^{XZs)d^ds 

Jo Jo 

-Ej 

r f G, (Xf A _ z,) d^ds. 

0 Jo 


In order to pass to the limit we shall need the following a-priori estimates. 
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First we apply the Ito formula for the functional u i— 5 > ^ \\u\^^ to the equation 


(4.7) 


dXl'^ = div dt + eXXl’^dt in (0, T) x O 

^ E B^{xt^)dt +j:b, (x!^^) dpi 

i—1 i—1 ^ ' 


+ - 


JCq — xq, 


dx: 


e,X 


= 0 , 


dv 

in order to get that 


in O 

on (0,T) X do 


(4.8) 


1 , 


-E 


X. 


£,A 


-E 


H 


[ [ {XXI’’')d^ds + eE [ [ d^ds 

Jo Jo Jo Jo 


< ^E\\xo\\l, VA, tG[0,r]. 


Moreover, in order to verify (13.31) . we see that by the Mosco convergence (see 
e.g. [33l Proposition 6.2]) 


= 


JO 


$ in Mosco sense as A ^ 0, 


and Fatou’s lemma (after passing to an a.e. convergent subsequence — strong 
L^-convergence is justified below), we get that 


E / $ 


0 


(X|) ds < liminf E [ [ <p^ d^ds 

'*'-*•0 Jo Jo 


< oo. 


On the other hand, also by the Ito formula and Lemma (i), we get that P-a.s., 

telo,T], 


^i^£jA2 


H 


+ 1^ (vn^i (VXf’^i) - 4/^^ (VXJ>^=), ds 


+ej^ 

= ^ / (B. - xj’^=) dl3l , xj’^^ - xj’^=)^ . 


By using [331 ec[- (A.. 6 ) in Appendix A], we have for all G and some 
positive constant C > 0 that 


- 'p^^(c),e - c) > -^(Ai + A2) (1 + lei" + Id"). 

We obtain P O ds-a.s. that 

(VX^’^i) - 4-^^ {XX^’^^) - VX®’^^)^ 2 (c,.Rd) 

> -C (Ai + A 2 ) y (1 + I VX^’^i 1^ + I VX^’^= 1^) d^ 

o 

> -C (Ai + A 2 ) (1 + |X"’^i |g + |X"’^= ly , 
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and then, by (|4.4I1 and (j4.8ll . we get for the expectation, that 

t 

E I ds > -C (Ai + A2). 

0 

Now, by the Burkholder-Davis-Gundy inequality, taking (12.21) into account, and 
by the above computation concerning 4''*', we obtain that for all t € [0,T], 

ifi sup + f IIVArj’^i - < ^(Ai + A2). 

2 0<s<i Jo 

Consequently, we have that 


(4.9) 

limE 

sup 

xD^ - xt 

2 


A-)-0 

iG[0,T] 


H 


= 0 . 


We can now pass to the limit for A —s- 0 in (14.61) in order to obtain (recall that 


(4.10) 


ifi \\X! -Zt\\l + E [' [ ^ {VXD d^ds 

+eE [ [ {VXD VXI - VZs) d^ds 
Jo Jo 

^Ellxo - Zoll^+E /" f (f{VZs)d^ds 
^ Jo Jo 

-E f [ Gs{Xl-ZDd^ds. 

Jo Jo 


Step IV (e —)• 0): 

Arguing as in the previous step, we get that 


(4.11) iE||Af||^+E r [ ^{VXDd^ds + eE f [ \VXI\^ d^ds 

^ Jo Jo Jo Jo 


< iEiixoll^, tG[o,r]. 


By Ito’s formula, this time, considering the process for fixed A > 0 and not £ > 0 
fixed, as previously, and by monotonicity, we get that 




-t 


H 


+ ^‘e-^ (£iVAJ^>^ - e2VAf VXf- VAf 

+ l f e-D\XD’^-XD’YHds 


, ds 


2 

N rt 


< E / dPl , AJ-^ - Af’^) 


H 
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Since P 0 ds-a.s., 

> -i (el IjAA'^i’^ll^ +ei ||AA^^’^||y - i 

and by using again the Burkholder-Davis-Gundy inequality and (12.21) . we get that 
for all t G [0,T], 

E sup e-'* - Af’^11^ 

0<s<i 

< efE [ e-" llAAJi’^ll^ds + e^E/" e”" || AAJ=’^||^ ds. 

Jo Jo 

Keeping in mind that for initial data in S), by Step II above, in particular, 

by ga, 

(4.12) eE f ||AAJ’^||^ds < C, 

Jo 

uniformly in A > 0, we obtain by (SH), for letting first A ^ 0, 


limE 

s —^0 


sup \\Xl 

te[o,T] 



= 0 


for some limiting process X G (7([0, T]; L^(n; id)). Note that by weak convergence 
in ij^([0, T];S), we get that P 0 dt-a.s. X G ii^(O). 

Finally, by computing P (8) ds-a.e. 

sf (VA^ - VZ) dC = -£ / AA^(A'^-Z)dC 
Jo Jo 

> -isl||AA^||^-ie§||A--Z||^ 

and using again (14.1211 we can pass to the limit in (14.1011 and get that A is a 
continuous SVI solution (which satisfies (j3.3p by passing to the limit in gUD) 


< 


\e\\x,-z,\\1 + e 


-E ||a;o — ^oll^ + IE 


f f (p (VXs) d^ds 

Jo Jo 

nt r 

p iXZs) d^ds 


0 JO 


-E 


0 JO 


Gs {Xs - Zs) d^ds. 


Step V (general initial conditions): 

In order to conclude the proof of existence we only need to extend the solution 
for arbitrary x G L^(f2, Aq,P; ii). Let A, A* be continuous SVI solutions starting 
in x,x*, resp. Note that S is dense in ii, so this follows directly from 


E 


sup ||At-A,*||^ 
te[o.r] 


< E ||a;o 


2 

H ’ 


which can easily be obtained by arguments similar to those from the previous steps 
(using monotonicity, Burkholder-Davis-Gundy and f|2.2p b 
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□ 


4.2. Uniqueness. 


Proof of Theorem \4- 1\ (continued). The existence of a continuous SVI solution is 

proved in the section above. We follow an argument from [33]. Let X be any 

continuous SVI solution to (El) with initial condition x G P; iL) and let 

ye,A,ra ^]^g gtj-Qjig approximating solution to (|4.5|) with initial condition y" G 

Lf{Pt,iFQ,T\S) such that y" —>■ y in where y G 

Then the following variational inequality holds (with Zq = y'^, Z = G = 

div 

^E\\Xt-Yf’^n\H+^l^ms)ds 

-E [ (eAy/’^’" + div4'^(VF/’^’”),A, 

Jo 

for a.e. t G [0, T]. 

By [33l Appendix A] for all z G we have 
-(div4'^(VF^’^’"),2-W’^’")ff+4>(y^’^’") < $( 0 )+C'A(l+$(y®>^’")) ds(g)P-a.e. 

Since $ is the lower-semicontinuous envelope of $ = (i.e., $ restricted to H^), 

for ds (g) P-a.e. (s,w) G [0,T] x fl, we can choose a sequence z™ G such that 
Xsioj) in H and $(z™) ^ 4>(As(a;)). 

Hence, 

-(div4'^(Vr^’^’"),X-W'^’”)ij+$(W'^’") < $(A)+C'A(1+$(F''’^’”)) ds®P-a.e. 
Thus, 

^EIIW - - yl\H + cxeJ^ (1 + $(y/'^’")) ds 

+ Ie^ (et ||Ay/’"’-|||, +£§||A, - ds. 

Taking first A —>■ 0 then e ^ 0 (using the bound (I4.12P which is uniform in A) 
and then n ^ oo yields 

EWXt-YtWlm^-yWl, 


for a.e. t G [0, T]. 


□ 


Appendix A. Vector fields of Killing 
In this section, suppose that O C open, bounded, with smooth boundary 

do. 

Definition A.l. A C^-vector field h : O ^ is called Killing vector field, if the 
following condition is satisfied on O 

(A.l) djP + diV =0 VI < i, j < d. 

The following lemma is based on ideas from [^, see also [55] . 
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Lemma A.2. A (sufficiently smooth) vector field b : O ^ is a Killing vector 
field if and only if 

(A.2) {Ab,b) = 0, div & = 0 onO. 

(A.3) ^ {djV + diK)vjb'^ = 0 on dO, 

where the Laplace acts componentwise and where v : dO —>■ denotes the outer 

unit normal of O. 

Proof. One easily sees that, 

dfiidjP + - b>{d,P)] 

={dAPW + (d.bWjb^) 

+ + {d,V){djP) 

-{d,V){d,K)-V{d,d,P). 

Assuming the above conditions, interchanging the order of differentiation, and 
summing over 1 < *, j < d, we obtain that 

^ + d^v)P -b>{d^P)] 

={Ab,b)+ [{d,Pr + id^lfi)im]-{dKb)^ 

l<i,j<d 

+ {djd^K)P - 

l<i,j<d 

E idjb^?+mbn{djb^) + {m^- 

l<i-,j<d 

= ^ E i^Jb^ + ^^Vy■ 

Denote by 5"^“^ the surface element on dO. By Gauss’s divergence theorem, we 
get that, 

I E [jdjb^ + df 

= E / [{djP+ d,V)P-V{d,P)]njdS‘^-^ 

= E / [{djP+ d^lf)P]iyjdS‘^-^ 

l<^,J<dd90 

= 0 , 

and hence djV + diV =0 on O. 

Suppose conversely, that 5 is a Killing vector field. Then (IA.3I) is automatically 
satisfied. Also, div6 = 0 by choosing i = j. Clearly, also 

0 = PdjdjP + Pdjd.V = b^djdjb^ + Pd^djbf 
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and summing over 1 <i,j < d yields, 

(A6,6) = 0, 

and hence (IA.2II is satisfied, too. □ 

Theorem A.3. Let b : O ^ be a -vector field. In order that the first order 
differential operator u i-A (6, Vu) commutes with the Laplace operator u i—> — Au on 
the space of smooth functions on O it is necessary and sufficient that b is a Killing 
vector field. 

Proof. See [33 Theorem 2.1]. □ 
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